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Abstract 

A percolation argument and a dilute compressible random field Ising 
model are used to present a simple model for mixed cyanide crystals. 
C^ ' The model reproduces quantitatively several features of the phase di- 

1^^ ■ agrams altough some crude approximations are made. In particular 

Q^ , critical thresholds Xc at which ferroelastic first order transitions disap- 

pear, are calculated. Moreover, transitions are found to remain first 
order down to Xc for all mixtures except for bromine, for which the 
transition becomes continuous. All the results are in full agreement 
'^ I with experimental data. 
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1 Introduction 

Mixed cyanide crystals X{CN)xYi^x where X is an alkali metal {K, Na or 
Rb) and Y stands for a spherical halogen ion {Br, CI or /) exhibit puzzling 
phase diagrams [1]. In particular the pure alkali-cyanide XCN ferroelastic 
transition disappears at some cyanide concentration threshold x^ which is not 
equal to the associated percolation threshold Pc- Moreover Xc varies with both 
X and Y compounds. For instance, Xc = 0.60 and Xc = 0.80 for respectively 

Y = Br and Y = CI with X = K. 

Below Xc experimental evidences hint to a cyanide orientational freezing, to- 
gether with no long range quadrupolar order. In parallel, the high temperature 
cubic symmetry is preserved at low temperatures [1] . Analogies to spin-glasses 
and random field systems were suggested to elucidate the orientational glass 
state physical nature [1, 2]. 

In addition, while the transition stays first order for X{CN)xCli-x [3], it 
was found to become continuous at Xc for X{CN)xBri^x [4]. These experi- 
mental facts are still lacking an explanation. 

The complicated symmetry of the molecules involved as well as the form of 
quadrupolar interaction make theoretical attemps rather difficult and heavy. 
Previous approaches [2 and references therein] tried to embody some of the 
molecule physical charateristics. However they failed to explain above features. 

In this letter a model which combines a dilute compressible random field 
Ising system with a percolation argument is suggested to reproduce parts of the 
mixed cyanide phase diagrams. A microscopic calculation of critical thresholds 
Xc is presented. Compressibility is found to produce a first order transition 
only at cyanide concentrations larger than a threshold x^. On the opposite, 
diluted random fields are shown to activate a first order transition solely below 
another cyanide threshold Xr- 

Transition orders which depend on the respective values of Xc, xl and Xr, 
are discussed for a large variety of mixtures. We predict transitions to remain 
first order down to Xc for all mixtures, except in the case of bromine which 
exhibits a continuous transition at Xc- Our results fit perfectly a series of 
experimental data. 

2 A percolation approach at T = 

Starting from the Xc 7^ Pc experimental fact, we conclude that some cyanides 
do not participate in the propagation of quadrupolar long range order. These 
cyanides are thus "neutralized" from reorientation. Due to steric hindrance, 
a local deformation of the unit cell is expected on substituting one spherical 

Y ion to the dumbbell-shaped cyanide. We assume that a local XY unit cell, 
embedded in XCN in bulk, matchs the unit cell of XY in bulk. The associated 
volume deformation of the former unit cell is Av = j{a^xcN ~ ^%x) ■ From 
symmetry, each one of the c Y nearest neighbors is affected by this local volume 



deformation which thus produces the overall volume deformation AV = cAv 
(c = 12 and pc = 0.198 on the fee cyanide sublattice [5]). 

AV^ may be extracted from neighboring unit cells using their free volumes 
which originate from volume differences between the unit cell and the molecule 
itself [6]. The free volume per molecule is Vf = -^f^ — {vx + vcn) , where 
o-xcN is the pure XCN lattice constant, vx and vcn are ion volumes [7]. 

Dumbbell-shaped CN reorientations are directly coupled to the unit cell 
shape via steric hindrance mechanisms. Any volume deformation of a cyanide 
cage (increase or decrease) lowers the corresponding symmetry which in turn 
lowers the number of accessible orientations. Here we are assuming that a 
cyanide whose cage is deformed (increased or decreased) becomes orientation- 
ally "neutralized". Accordingly, one Y substitution will affect on average, 

a = c , (1) 

unit cells by deforming their respective free volumes. As a consequence, we 
obtain an effective density of free to reorient cyanides, Xf = x — a{l — x) . 
These free to reorient cyanides thus obey, by definition, site percolation with 
Xf^c = Pc at X = Xc- Therefore we get, 

.0 = ^. (2) 

1 + a 

Without a fitting parameter the calculation of Xc is readily performed using 
crystallographic data [7, 8, 9]. The results are obtained for various mixtures 
(see the Table) including systems for which no experimental data are avail- 
able. Below Xc a region of randomly oriented ferroelastic domains with no 
static phase transition is predicted. These domains will shrink with increas- 
ing dilution to disappear eventually when x/ = at a new threshold [10], 
Xd = YT^ • Shear torque experiments [11], as well as diffraction experiments 
[12] suggested a very similar phase diagram. 

3 A crude Ising-like model 

To extend above results to T 7^ we now build the simplest possible model 
Hamiltonian making several crude approximations. 

• From symmetry, cyanides have several equivalent orientations making 
g-Potts variables appropriate [13]. However we restrict ourselves to the 
minimum number of orientations required to sustain an orientational 
long range order, i.e., two. Therefore we use Ising variables {Si = ±1} 
to mimic the cyanide orientational degrees of freedom. 

• Cyanide quadrupolar interactions are long ranged. However our per- 
colation argument is a short range effect. Therefore, to be consistent 
with above T = calculation, we consider short range ferromagnetic 
couplings. 



To fit the first order character of the pure XCN transition, we introduce 
elastic degrees of freedom in order to turn the continuous Ising ferro- 
magnetic transition to first order. To keep calculations simple, we use 
an harmonic model of volume fluctuations though it is clearly a rudi- 
mentary model of elasticity [14]. 

Within a model of two equivalent orientations, a cyanide prevented from 
reorientation is trapped along one direction. To embody this effect, we 
introduce local quenched random fields. The probability pt to have a local 
random field is equal to the probability of having a deformed cyanide 
cage. From the density Xt of trapped cyanide we obtain, 

( a(l-x) ]f ^^ < ^ < 1 

Pt=\ ,- ! ' - - • (3) 

[1 lix < Xd 

The distribution function for the random filed hi is then, 

P{hi) = ^[5{h, -h) + 6{h, + h)] + (1 - pt)6{h,) . (4) 

It satisfies both required symmetry conditions P{hi) = P{—hi) and hi = 
(the overline denotes a field configurational average) which preserve 
the cubic symmetry of the associated problem. 



4 The Hamiltonian 

From above approximations, we obtain the following effective Hamiltonian, 

E 
Heff = -G J2 ^i^jSiSj - ^( XI ^i^^jSiSjf - Y^ eihiSi , (5) 

where G and E are constants, A^ is the total number of spins, d is the dimension 
[15], and e^ is a random site variable. It is 1 if site i is occupied by a cyanide 
and otherwise. We have {ej}a„ = x, where {...}av denotes a configurational 
average over site disorder. 

We define the order parameter as m = {< tiSi >}av where averages are 
taken over both site and random field disorders. The associated mean field 
site free energy is, 

JF = -cGm^ H m^ 

2 4 



— {ln[cosh{(3cGm + [3Lnt' + (3h)] + ln[cosh{(3cGm + [3Lm^ - (3h)]} 



—xksT 

L z 

-f(l - pt)ln[cosh{f3cGm + /SLm^)] 



kBTln{2) , (6) 



where (3 = ^^, ks is the Boltzman constant, T is the temperature and L = 
c^E (more details will be published elsewhere). In real systems hi and Cj are 



correlated. However, in order to be consistent with the previous short range 
interaction approximation, to keep the calculations simple and doing a mean 
field calculation, these correlations are neglected here. We now analyse two 
simple limiting cases which are physically meaningfull. 

4.1 The zero-steric hindrance effect case {a = 0) 

From a Landau expansion of Eq. (6) a continuous transition occurs at the 
critical temperature, ksTc = xcG under the condition of a positive quartic 
coefficient, B = —-^ + -^ , which results in the condition x < xl on cyanide 
concentration [16] where, 

^L - i^Y" . (7) 

Ax, X = Xl the transition turns first order via a tricritical point (5 = with a 
positive free energy sixth order coefficient). 

At X = 1, Pi = (Eq. (3)) even if a 7^ 0. From experimental works 
pure XCN exhibits a first order transition. All plastic systems must therefore 
satisfy x^ < 1 with thus L > ^. Dilution weakens the first order character 
of the transition. The associated negative quartic term B becomes smaller 
in amplitude to vanish eventually at a tricritical point {x = xl)- There, the 
transition is continuous with tricritical exponents. Upon further dilution the 
transition becomes second-order with < x < X/,. 

4.2 The zero-compressibihty case (L = 0) 

In this case, random field are acting alone (a 7^ 0). At x = 1, with pt being an 
independent external parameter, the zero-compressibility free energy becomes 
identical to that of the trimodal random field Ising model [17]. A first order 
transition is found only for 0.73 < pt < 1 and for some restricted range of 
random field intensities ~ 0.55 < -^ <~ 0.65 [17]. These results produce an 
additional threshold in cyanide density around pt = 0.73 (Eq. (3)), 



0.73 + a 

Only dX X < Xj. can dilution turn the transition to first order via random 
fields. However condition Xc < x^ must also be satisfied since the transition 
itself disappears at Xc. The equivalent constraint on a gives, 

OTSpe 2.09 . (9) 

1-0.73-pc 

Eq. (8) gives Xr = 0.60 and x^ = 0.83 for respectively K{CN)xBri-x 
{a = 1.12) and K{CN)xCli^x (« = 3.48). In parallel the Table shows Xc = 
0.62 and Xc = 0.82 for Br and CI respectively. On this basis we conclude 
that upon dilution, while random fields can turn the transition first order 
in K{CN)xCli-x mixtures (xc < x^), they cannot do it for K{CN)xBri^x 
mixtures (xc > x^). 



5 Conclusion 

In our model, both compressibility and random field drive the transition to 
first order. However, the compressibility effect weakens with dilution to fade 
out in the vicinity of Xc for all systems. On the opposite, random fields start 
to be active at dilution below Xr, which is possible only when Xr > Xc- Using 
the condition a > 2.09 (Eq. (9)), the Table shows it is always satisfied except 
for mixtures with bromine. 

We can thus predict that dilution with chlorine and iodine maintains the 
transition first order down to Xc- Only bromine turns the transition continu- 
ous. A tricritical point is expected for X{CN)xBri-x mixtures. Our predic- 
tions reproduce experimental results with respect to potassium systems [3, 4]. 
Additional experiments on mixtures with sodium and rubidium would provide 
a definite ground to our model. 

It was worth to stress that the extension of the model to g-Potts variables 
will not discard elasticity with respect to the first order character of the tran- 
sition. Indeed, the rigid Potts model does exhibit a first order transition at 
d = 3 , but only when q > 2. It is the case for KCN, but not for all sys- 
tems. In particular, most materials which have a tetrahedral molecule in a 
cubic site will have equivalent orientations with Td symmetry, i.e., q = 2, like 
for instance the perchlorate tetrahedron in KCIO^. Therefore, elasticity is a 
necessary feature of our model. 

Last but not least, the predictive power of our model comes as a surprise 
according to the several approximations made. In particular the short range 
character of the interactions is hard to justify. However, the results are impres- 
sive showing there must exist some screening mechanism which allows such a 
simple model. 
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Table 1: Numerical values calculated for a, Xc and Xd and experimental thresholds 
when known (denoted by 'exp:'). Units for lengths and volumes are A and A^. 
See details in the text. Error bars are within respectively ±0.01 for all data in the 
Table, and ±0.05 for experimental thresholds. 



XCN/Y 


axcN 


axY 


Av 


Vf 


a 


Xc 


Xd 


KCN/Cl 


6.53 


6.29 


7.30 


25.19 


3.48 


0.82 exp: 0.80 


0.78 exp: 


0.75 


KCN/Br 


6.57 


6.60 


2.36 


25.19 


1.12 


0.62 exp: 0.60 


0.53 exp: 


0.50 


KCN/KI 


6.53 


7.06 


18.46 


25.19 


8.79 


0.92 exp: 0.90 


0.90 




NaCN/Cl 


5.90 


5.65 


5.99 


12.67 


5.68 


0.88 exp:r^ 0.80 


0.86 




NaCN/Br 


5.90 


5.97 


2.11 


12.67 


1.99 


0.73 


0.67 




NaCN/KI 


5.90 


6.47 


16.63 


12.67 


15.74 


0.95 


0.94 




RhCN/Cl 


6.82 


6.58 


8.05 


30.98 


3.12 


0.80 


0.76 




RhCN/Br 


6.82 


6.85 


1.19 


30.98 


0.46 


0.45 exp: 0.55 


0.32 





